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Abstract 

The  reduced  size  and  operational  energy  associated  with 
microelectronic  circuitry  has  created  a  situation  wherein  a 
cosmic  ray,  in  a  single  pass,  can  generate  a  false  signal  or 
otherwise  damage  a  circuit  element.  In  this  work  an  analytic 
expression  for  the  event  rate  of  false  signals  (a.k.a.  soft 
errors)  is  derived,  based  upon  the  geometric  probability  of  a 
cosmic-ray  track  occurring  in  the  sensitive  volume  of  a  micro- 
circuit  device.  The  analyses  show  that  substantial  reduction 
in  the  soft  error  rates  from  heavy  cosmic  rays  can  be  achieved 
by  device  designs  which  incorporate  large  operational  energies 
and  small  depletion  collection  volumes.  The  permanent  damage 
to  oxide  layers  in  metal -oxide-semiconductor  (MOS)  devices  is 
estimated.  The  special  role  of  protons  as  initiators  of  soft 
errors  and  permanent  damage  is  explored.  These  processes  are 
then  described  in  the  context  of  very  large  scale  integration 
(VLSI)  electronics  in  Earth  satellite  environment.  Proton 
non-elastic  scatter  events  which  produce  energetic  recoiling 
Silicon  nuclei  are  identified  as  the  chief  cause  of  proton 
initiated  soft  errors. 


Introduction 

Microelectronic  circuits  (VLSI)  subjected  to  cosmic-ray 
or  other  heavy  ion  bombardment  undergo  false  switching-error 
writing  response.  The  requirement  of  the  circuit  to  produce 
such  a  response  is  the  deposition  of  a  threshold  amount  of 
energy  in  the  form  of  ionizing  reactions  in  a  small  fairly 
well-defined  volume  of  the  circuit.  The  requirements  for  the 
ion  or  cosmic  ray  are  that  the  energy  be  sufficient  to  pene¬ 
trate  any  shielding,  with  kinetic  energy  in  excess  of  the 
circuit  threshold,  AE,  and  linear  energy  transfer  great  enough 
to  deposit  AE  in  a  track  whose  length  is  only  several  microns. 

This  work  describes  the  dose  distribution  surrounding  a 
cosmic- ray  track  and  the  impact  of  that  energy  deposition  in 

This  paper  is  declared  a  work  of  the  U.S,  Government  and 
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the  form  of  hole-electron  pairs  in  VLSI  circuitry,  specifically 
as  it  applies  to  dynamic  Random  Access  Memories.  An  event 
rate  for  errors  in  satellite  borne  electronics  exposed  to  the 
cosmic  flux  is  derived.  The  possibilities  of  permanent  damage 
to  the  oxide  layers  in  MOS  devices  are  described,  and  finally 
the  role  of  protons  in  these  interactions  is  developed. 

Size  and  Dose  Distributions  near  Heavy  Ion  Tracks 

When  a  fast  heavy  ion  penetrates  condensed  matter,  some 
of  the  orbital  electrons  are  stripped  from  the  ion,  and  a  beam 
of  such  ions  then  has  a  mean  ionic  charge,  z  <_  Z;  z  has  been 
measured  widely  as  a  function  of  energy  and  velocity  of  the 
ion  and  is  given  by  Barkas1 

z  ■  Z(1  -  e-lse2'2/3>  (1) 


The  ion  loses  energy  through  electron  collisions  and  generates 
electron  flux  at  a  rate  given  by  Mott2 


dn  2ttNz26^ 

doj  2 

mv 


i  -  »2  “'n»  *  Br 


(2) 


where  dn/dw  =  the  number  of  electrons  between  energy  oi,  ai  +  doj 
generated  per  cm  of  ion  track  In  matter  containing  N  electrons/ 
cm3  and  tom  =  2mc2$2Y2;  e,  m  are  the  electron  charge  and  mass 
and  8  B  v/c. 

The  above  electron  generation  mechanisms  can  be  coupled 
with  electron  transport  models  to  produce  a  relatively  detail¬ 
ed  picture  of  the  dose  and  electron. flux  near  a  heavy  ion 
track.  Kobetich  and  Katz3  and  Hamnr  have  made  extensive  cal¬ 
culations  of  the  energy  deposited  near  an  ion  track  with  the 
following  result: 

D  *  kz2/r2  (3) 

That  is,  the  energy  deposition,  D,  In  ergs/g  at  radius  r  from 
a  cylindrical  track  falls  off  as  1/r2  and  is  proportional  to 
z2,  where,  from  Eq.  1  above,  z  is  the  mean  Ionic  charge  of 
the  ion  or  cosmic  ray.  Both  z  and  k  are  functions  of  8.  Re¬ 
cent  measurements  at  Brookhaven®  have  verified  the  1/r2  depend¬ 
ence  over  a  range  of  1  to  300  nm  (10  A),  as  shown  In  Fig.  1. 

The  column  of  ionization  which  a  fast  heavy  ion  leaves 
in  its  wake  contains  central  core  energy  densities  of  8xl020 
eV/cm3  for  protons  to  103  times  greater  for  heavy  cosmic  rays 
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(e.g.,  Fe88).  The  track  can  be  viewed  as  a  cylinder  whose 
radius  is  approximately  lu  and  within  which  the  dose  falls 
off  from  the  center  as  r"‘.  For  this  work  it  is  the  size  of 
the  cylinder  in  relation  to  the  size  of  electronic  devices,  as 
shown  in  Fig.  2,  and  the  energy  deposited  therein  that  is  of 
interest. 


Track  Length  Distribution  in  3-Dimensions 

Recent  increases  in  the  density  of  planar  technology 
electronic  microcircuits  have  led  to  devices  which  are  capable 
of  interruption  by  the  passage  of  a  single  cosmic  ray  or  other 
highly  ionizing  particle.  One  of  the  more  prominent  of  these 
is  the  dynamic  random  access  memory  cell  of  the  MOS  type.  The 
interaction  of  fast  heavy  ions  with  MOS  memory  cells  can  create 
writing  errors.  Two  sources  have  been  identified:  1)  cosmic 
rays6"8  and  2)  uranium  and/or  thorium  a  decay9  from  ceramic 
packaging  materials  adjacent  to  the  memory  cell.  The  under¬ 
lying  requirement  for  memory  cell  upset  Is  the  deposition  of 
an  amount  of  energy,  greater  than  some  threshold  energy,  which 


Fig.  1  Comparison  of  radial  distributions  for  heavy  Ions, 
SEM  electron  beams,  theory,  and  experiment. 
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Fig.  2  Heavy  ion  track  in  perspective  of  an  MOS  device  of 
VLSI  size.  Concentric  circles  are  dose  in  rads  in  silicon. 


generates  an  electron-hole  pair  flux.  The  subsequent  collect¬ 
ion  of  the  charges  causes  a  writing  error.  The  energy  must 
be  deposited  in  the  sensitive  region  of  the  device,  i.e., 
where  the  charges  are  collectable.  In  general  that  region  is 
near  the  gate  and  source  of  the  MOS  device.  The  possibility 
of  large  diffusion  lengths  enlarges  that  region.  The  dose 
surrounding  a  cosmic-ray  track  as  it  impinges  normally  on  an 
MOS  structure  is  shown  in  Fig.  3. 

Analysis 

The  approach  to  the  problem  taken  here  is  to  treat  the 
collection  volume  (sensitive  volume)  of  the  electronic  device 
as  a  right  rectangular  volume.  The  cosmic  rays  traverse  this 
volume  without  deflection,  and  the  tracks  are  thus  straight 
lines.  This  permits  a  solution  based  on  the  geometric  proba¬ 
bility  of  the  chord  lengths  generated  when  a  convex  volume 
(surface  everywhere  characterized  by  a  positive  radius  of 
curvature)  is  immersed  in  a  uniform  isotropic  flux. 

An  ion  moving  through  matter  loses  energy  at  a  rate 
called  its  LET  (linear  energy  transfer).  For  small  distances 
this  loss  rate  Is  a  constant.  The  total  energy  lost  by  the 
Ion  is  the  LET  x  l,  where  I  Is  the  track  length  in  the  region 
of  Interest.  The  assumptions  of  this  work  are  that  the  heavy 
ion  tracks  In  question  are  straight  lines  and  that  the  volume 
of  Interest  is  a  rectangular  parallel  piped,  hxaxb,  h<a 
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SOURCE  GATE  DRAIN 

Fig.  3  Dose  distribution  surrounding  a  heavy  ion  track 
normal  to  VLSI  device. 


<  b.  Actual  typical  dimensions  are  h*2p;  a,b -  20p  in  large 
scale  Integration  (LSI)  and  h,  a,  b  <  lp  in  VLSI. 

What  is  required  Is  the  probability  that  a  cosmic-ray 
track  length,  within  the  collection  volume,  is  greater  than 
or  equal  to  the  minimum  length  needed  for  the  electronic  re- 
sponse  to  occur.  The  probability  distribution  which  answers 
that  requirement  is  called  the  sum  distribution.  The  sum 
distribution  for  chord  lengths  in  a  rectangular  volume,  C(£), 
was  determined  exactly  from  the  work  of  Coleman10  who  solved 
the  problem  for  the  2-dimenslonal  rectangle  and  unit  cube,  of 
Kellerer11  who  developed  the  formalism  which  leads  to  the  3- 
dlmenslonal  sum  distribution  for  right  cylinders  of  arbitrary 
cross  section,  and  of  Bradford12  who  applied  the  formalism 
to  a  rectangular  parallel  piped. 

The  resulting  analytic  functions  C(t)  are  plotted  In 
Figs.  4  and  5  for  several  cases.  The  symmetry  of  the  infinite 
plane  slab  and  the  sphere  permits  easier  evaluation  of  the 
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integrals,  and  they  are  shown  for  references.  The  distrib¬ 
ution  for  the  unit  cube  is  shown,  as  is  the  distribution  for  a 
volume  h  x  4h  x  6h.  The  latter  curve  is  pertinent  to  the 
attempts  of  Pickel  and  Blandford8  to  predict  an  event  rate 
for  cosmic-ray  interruption  of  satellite  borne  electronics. 

The  volume  in  question  is  the  sensitive  volume  of  each  of  24 
4K  RAM  memory  cells,  generally  that  volume  defined  by  the 
source- to-drain  area  and  depletion  depth. 

We  can  compute  an  event  rate  for  electronic  interrupt 
using  C(i.)  by  noting  that  any  convex  body  immersed  in  a  uni¬ 
form,  isotropic  flux,  9  parti  cl es/cmz/sec,  will  experience 
(«pS/4)  transits/s,  where  S  is  the  total  surface  area. 

♦(x),  the  LET  spectrum,  specifies  the  particle  flux  for 
each  value  of  LET.  A  minimum  track  length  along  which  suf¬ 
ficient  energy  is  deposited  to  cause  the  electronic  response 
to  occur  is  determined  by  each  LET  value.  The  relationship 
which  leads  to  that  minimum  track  length  is 

l  =  ^threshold 
LET-p 

where  p  is  the  density  and  the  AE  threshold  is  the  minimum 
energy  deposition  required  for  electronic  response.  C (£) 
then  gives  the  fraction  of  track  lengths  generated  by  the 
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Fig.  5  Differential  and  sum  distributions  for  track  lengths 
in  VLSI  device.  Volume  is  of  rectangular  shape,  h  x  4h  x  6h. 

particle  flux  with  that  LET  value  which  can  cause  the  elec¬ 
tronic  response.  Integration  over  the  LET  spectrum  then 
yields  the  total  event  producing  combinations  of  LET  and 
sufficient  track  length.  The  event  rate  thus  is  given  by 

n  s  S/4  / *1  <t(x)  C(AE/px)dx  events  per  sec  per  device  (4) 


assuming  100%  efficiency  for  the  process  (e.g.,  charge  collect¬ 
ion),  4(x)  is  the  LET  flux  in  particles/fcm?  .s-(MeV  cmvg)] , 
x  is  LET  in  MeV  cm2/g  and  <p  =  /$(x)dx 

Since  C (£.)  is  a  reasonably  complicated  function,  the 
event  rate  integral  usually  requires  machine  solution.  An 
approximation  can  be  obtained  by  noting  that  the  dominant 
terms  in  C(&)  fall  off  like  l/ir,  r  >  2  and  that  4»(x)  -  4>q 
(xo/x)p  •  P  =  3.8. *3  Equation  (4)  is  thus  reducible  to  simple 
form  by  taking  the  following  definitions  for  the  integrand: 

*(xHx0)(x0/x)p=i6xio13(io3/x)3.8 
part1c1es/[m2«  s  •  (MeV  cm2/g)]  (5a) 

with  the  values  for  x0,  4>(xo),  P  taken  from  Ref.  13.  For  C(t) 
take 

CU)  =  C(h)(h/*)2*2  «  0.75  (hpx/AE)2,2  (5b) 
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as  a  reasonable  approximation  for  chord  lengths  greater  than  h 
in  rectangular  volumes.12 

x0  =  AE/p)W  =  1q3  MeV  cn,2/9  (6a) 

Xj  =  AE/pS-m^n  =  7xl08  Mev  cm2/g  (6b) 

are  the  limits  in  the  integral  of  Eq.  (4),  where  the  numbers 
are  for  the  4K  RAM  of  Ref.  8.  The  effect  of  these  limits  in 
determining  the  portion  of  two  species  C  and  Fe  of  the  cosmic 
spectrum,  which  can  participate  in  that  case,  is  shown  in 
Figs.  6  and  7. 

The  general  solution  is  obtained  by  retaining  the  defi¬ 
nitions  of  the  limits  following  integration  of  Eq.  4  and  the 
event  rate,  n,  is  given  by 

n  =  7.5xl0‘7NS(hp/AE)2'2l(x0)'-6-(x1)”*6] 

=  7 . 5x10'  7NS  ( hp/ AE ) 2  • 2  [( AE )  • 6-  ( ptmi  n/ AE )  * 6  ] 

=  7.5xlO"7NS(hp/AE)2"8[(fc|nax/h)"6  -ll 
»  7.5xl0'7NS(hp/AE)2,8[(10AE/hp)*6-l]  events/ 


I  10  10*  10*  To4 

MeV  / NUCLEON 


Fig.  6  Cosmic  flux  for  Fe56  with  LET  values  for  the  same 
energy  range.  Shaded  area  shows  population  which  can  cause 
electronic  Interrupt  of  Ref.  8. 
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MeV  /  NUCLEON 

Fig.  7  Cosmic  flux  for  with  LET  values  for  the  same 
energy  range.  Shaded  area  shows  population  which  can  cause 
electronic  Interrupt  of  Ref.  8. 

Since  &min  *  h  in  this  approximation,  Umax  is  the  major  diago¬ 
nal  of  the  volume.  This  result  is  valid  for  cosmic  ray  flux 
with  LET  >  103  MeV  cmvg,  where 

.  .  2 

S  =  surface  area  of  device  sensitive  volume  in  u 

h  =  device  sensitive  volume  minimum  dimension  (presum¬ 
ably  depletion  depth)  in  u 

AE  =  switching  energy  in  MeV 
N  =  number  of  devices 
p  =  density  in  g/cm3 

This  result  is  based  on  the  unshielded  LET  spectrum  of 
Heinrich*3  for  6<z<26.  To  take  shielding  into  account  one 
need  only  multiply  these  results  by  an  attenuation  factor 
from  Heinrich's  paper,  -e*3,91(t/50) ,  where  t  is  shield  thick¬ 
ness  in  g/cm^  This  factor  applies  only  to  the  total  LET 
spectrum  >  103,  not  to  any  individual  species. 

The  event  rate,  Eq.  (7),  can  be  used  by  electronics 
designers  or  system  analysts  to  anticipate  soft  error  rates  in 
satellite- borne  electronics  devices.  The  basic  requirement 
for  such  utilization  is  to  identify  those  memory  cells  or 
switches  which  can  be  characterized  as  having  a  rectangular 
collection  volume  and  within  which  a  threshold  energy  deposi¬ 
tion,  AE,  can  be  specified.  By  way  of  example,  when  the 
numbers  from  the  Pickel  and  Blandford^  analysis  of  the  4K 
RAMs  in  satellites  are  used  in  Eq.  (11),  one  gets  203  events/ 
year  for  the  98,309  cells  involved.  This  compares  reasonably 
well  with  the  rate  calculated  in  Ref.  8. 
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Equation  (7)  has  further  merit  in  that  it  demonstrates 
the  event  rate  sensitivity  to  physical  parameters.  Clearly 
emphasized  is  the  necessity  to  scale  the  minimum  dimension,  h 
(depletion  depth),  faster  than  the  threshold  energy,  &E,  as 
one  designs  the  higher  density  VLSI  technology.  Failure  to 
utilize  an  analysis  which  includes  the  track  length  distrib¬ 
ution  can  be  seen  in  Ref.  8  where  the  mean  chord  length  is 
used  instead  of  the  full  distribution.  This  leads  to  the 
erroneous  conclusion  in  that  paper  that  the  minimum  LET  rate 
which  can  produce  the  electronic  error  is  4.6x103  MeV  cm2/g, 
whereas,  in  fact,  it  is  9.6x10?  MeV  cm?/g.  Finally,  as  is 
intuitively  clear,  the  longer  track  lengths  can  come  only 
from  ions  with  larger  angles  of  incidence  to  one  of  the  planes 
of  the  device.  This  response  to  angle  of  incidence,  which 
agrees  with  the  behavior  of  C(£),  has  been  measured  exneri- 
mentally.l4’!5  Thus,  establishing  a  threshold  energy  for  the 
device  error  sets  a  limit  on  the  angle  of  incidence  for  each 
LET  value. 


Permanent  Damage  (Oxide) 

Permanent  damage  in  the  oxides  of  MOS  devices  occurs 
when  charges  are  trapped  in  the  oxide.  That  charge  presence 
causes  a  shift  in  the  threshold  voltage  and,  with  sufficient 
shift,  a  disabling  of  the  device.  Normally,  this  effect  is  of 
concern  to  nuclear  survivability  and  is  caused  by  the  y  and 
x-ray  fluxes  of  that  environment.  Nonetheless,  a  cosmic  ray  in 
a  single  pass  can  create  enough  charge  pairs,  which  are  sub¬ 
sequently  trapped,  to  cause  threshold  shifts.  In  VLSI  the 
voltage  shifts  required  to  disable  the  device  are  smaller,  and 
hence  one  can  anticipate  that  a  permanent  damage  rate  will 
ensue. 

Equation  (7)  can  be  used  to  some  value  in  this  case.  AE 
is  set  as  the  energy  required  to  achieve  some  average  dose  in 
the  oxide;  then 

AE(MeV)  =  1.6xl08D(rads)V(cm3)  (8) 

If  the  area  is  shrunk  to  VLSI  dimension  with  V  =  O.lxlxlxlO*^ 
cm3,  for  104  rads  or  greater,  AE  =  0.16  MeV  and  n  =  1  events/ 
year/lO^  devices;  for  103  rads  or  greater,  AE  =  .016  MeV,  but 
EMminP  =  606  MeV  cm?/g,  which  violates  the  validity  criterion 
(LET  >_  103)  for  Eq.  (7).  Therefore,  one  must  resort  to  the 
formal  solution  of  Eq.  (4)  over  the  entire  LET  range  for 
evaluation.  Many  times  as  many  ions  qualify  as  a  result  of 
this  lower  LET  requirement. 

One  should  keep  in  mind  that,  because  of  the  dramatic 
nonuniformity  of  the  dose  distribution  in  the  track,  even 
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though  the  average  dose  is  greater  than  10^  rads,  significant 
regions  of  the  oxide  have  received  >  105  rads.  The  response 
of  devices  to  such  nonuniform  irradiation  is  not  well  es¬ 
tablished. 

Expected  Contributions  to  LSI/VLSI  Soft  Error 
Due  to  Solar  and  Cosmic  Protons 

The  interaction  of  heavy  ions  with  microelectronics 
devices  was  examined  in  the  preceding  section.  The  nature  of 
the  interaction  is  such  as  to  produce  error  signals  or  switch¬ 
ing,  and  the  chief  means  of  avoiding  those  errors  seems  to 
be  1)  one  very  small  dimension  associated  with  the  sensitive 
volume  of  the  device,  2)  large  AE  for  writing  or  switching, 
and  3)  redundant  programming  or  circuits.  In  all  the  fore¬ 
going  analyses  the  mode  of  energy  deposition  is  the  LET  due 
to  ionizing  collisions.  The  foregoing  analyses  also  recog¬ 
nize  a  minimum  in  the  LET  rate  necessary  to  produce  the  un¬ 
desirable  effects.  That  minimum  is  currently  about  1(H  MeV 
cm2/g,  and  the  cosmic  flux  greater  than  that  value  has  been 
characterized^  by  a  simple  expotential  [e.g.,  d>jf LET )  = 
fcoUo/xr  •  *0  =  4xl0"3m_2  srd"l  sec"l  MeV-1  cm-2  g}. 

The  vast  majority  of  the  cosmic-ray  population  has  LET 
rates  less  than  10' 3  MeV  cm2/g  however.  Protons,  for  example, 
have  maximum  rates  of  only  '102  MeV  cm2/g  at  energies  of  a 
few  MeV.  Since  the  preponderant  portion  of  the  proton  cosmic 
flux  has  LET  rates  of  2  MeV  cm2/g  (480  ev/p  in  Si)  for  a  lp 
size  device,  the  threshold  energy,  (AEth),  would  need  to  be 
as  small  as  500  eV  for  protons  to  qualify  for  error  production 
via  ionization  loss.  Protons  can  participate  even  when  AEth 
is  large,  however,  via  nuclear  reactions  and  by  elastic 
scattering. 

Thus,  the  role  protons  play  must  be  examined  in  two 
views:  1)  in  trapped  belt  or  solar  flare  exposures,  in  spite 
of  the  small  LET  rates,  the  flux  can  be  large  enough  to  cause 
permanent  oxide  damage  through  accumulated  total  dose;  2) 
cosmic  and/or  trapped  belt  plus  solar  flare  exposure  with 
protons  acting  singly  to  produce  soft  errors.  The  permanent 
oxide  damage  is  not  different  from  that  encountered  in  x-ray 
fluxes  from  nuclear  weapons  and,  since  studies  of  that  effect 
are  plentiful  in  literature,  it  will  not  be  dealt  with  here. 

Nuclear  Interactions 

The  proton- induced  nuclear  reactions  of  interest  are  p, 
n;  p,  2n;  p,  pn;  p,  ot,;  and  p,  spall.  The  mechanism  by  which 
these  reactions  can  trigger  soft  errors  is  the  ionization 
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energy  loss  of  the  recoil  residual  nucleus  and,  in  the  case  of 
the  p,  a,  and  p,  spall  reactions,  through  the  ionization 
losses  of  the  reaction  products.  The  total  proton  inelastic 
scatter  cross  section  rises  from  a  few  hundred  millibarns  at 
4  MeV  to  1200  mb  at  30  MeV  and  then  declines  to  -100  mb  over 
several  hundred  MeV,  averaging  350  mb  from  100  to  400  MeV.’8*!7 
The  p,  p';  p,  n;  p,  2n;  p,  n  reactions  contain  the  majority  of 
this  cross  section.  For  these  reactions  to  trigger  soft 
error,  the  recoiling  target  nucleus  must  be  the  agent  for 
ionization  energy  deposition.  Unfortunately,  the  recoil 
kinetic  energy  distribution  has  not  been  information  of  acute 
interest  to  nuclear  researchers  and  often  is  lumped  into  the 
residual  excitation  energy.  Thus,  the  fraction  of  the  in¬ 
elastic  reactions  which  lead  to  a  recoiling  nucleus  of  suf¬ 
ficient  energy  to  create  soft  error  is  not  readily  available. 

Low-energy  a  particles  already  have  been  identified  as  a 
cause  of  soft  errors  in  RAM  cells.  Thus  one  is  bound  to 
investigate  reactions  which  produce  a  or  larger  fragments. 

The  Si  (p,  a)  reactions,  as  with  all  light  nuclei,  have 
relatively  flat  cross  sections  for  Ep  100  MeV  at  about  12  mb 
after  peaking  at  60  mb  at  25  MeV.  ’  18  The  Q  value  for  He8 
production  is  16.9  MeV  and  for  He^  is  7.7  MeV.  At  proton 
energies  several  times  these,  the  p,  a  is  a  direct  reaction, 
e.g.,  little  residual  nucleus  excitation  or  recoil  and,  as  a 
result,  the  a  particles  are  produced  in  the  forward  direction 
at  high  energies  (Ep  -  Q).  In  this  energy  range  particles 
will  have  LET  rates  of  the  order  of  5  keV/p  and,  in  the 
dimensions  of  LSI  (20p),  cannot  deposit  the  required  energy 
(-5  MeV).  It  would  require  moderation  (i.e.,  through  shield¬ 
ing)  of  the  incident  proton  flux  to  energies  in  the  tens  of 
MeV  range  to  produce  an  abundance  of  a  particles  with  low 
energies  and  high  LET  rates. 

Spallation  reactions  have  been  measured^  in  Al^7  and 
reflect  total  cross  sections  of  about  6  mb  for  fragment 

production  (i.e.,  F^,  o^9,  Be7).  Although  the  fragments  may 
differ  with  a  Si 28  target,  the  strength  of  the  cross  section 
will  not  change  much  from  Al?7.  Thus,  the  combined  p,  a  and 
p  spall  cross  section  ( -  17  mb  for  Ep  >  100  MeV)  leads  to  the 
following  estimated  reaction  rates: 

n  reactions/year  =  Voton/yr*cm2xpNAV  {9) 

In  the  LSI  case  analyzed  by  Pickel  and  Blandford8  V  =  1030p3 
and  n  »  2.3xl0-fy(year*device)  =  23/(year*105  device)  due  to 
p,  a  and  p  spall . 
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In  a  solar  flare,  proton  flux  levels  increase  dramatic¬ 
ally  by  factors  of  as  high  as  10$  for  periods  of  the  order  of 
24  hr.  In  such  a  flare  (10^,  24-hr  duration)  the  total  num¬ 
ber  of  p,  a  and  p  spall  reactions  in  a  10$  RAM  would  be  -6300. 
By  using  the  total  average  inelastic  cross  section  (350  mb), 
one  would  see  460  inelastic  scatters  per  year  per  105  devices 
for  the  cosmic-ray  flux  and  140,000  reactions  in  a  solar 
flare  of  10^  x  cosmic  flux  lasting  for  24  hr.  The  number  of 
these  which  lead  to  soft  errors  is  not  easily  known  and  re¬ 
quires  the  energy  spectrum  of  the  reaction  products  for 
solution.  Although  the  event  rate  is  small  for  each  VLSI 
element  (since  the  volume  is  - 10"3  that  of  LSI),  this  is  gen¬ 
erally  compensated  for  by  the  increased  number  of  VLSI  ele¬ 
ments  per  chip. 

Elastic  Scattering 

Consider  a  proton  of  kinetic  energy  100  MeV  incident  on  a 
Si  nucleus,  mass  number  28.  Assume  that  a  transfer  of 
greater  than  100  keV  is  necessary  by  elastic  collision  for 
soft  error  to  occur.  The  cross  section  is  obtained  from  the 
integrated  Mott-Rutherford  relation  for  energy  transfer 
greater  than  Ep. 

ag  =  ( TTb^/4y^ ) { ( e- 1 )  -  B2lne+iraB  [2(e^-l)  -  lne]}  (10) 

where 

y  =  (1  -  B2)'?s  B  =  v/c 

c  -  VE0  «  '  22e2/* 

E  =  2£(E  +  2mc2)  . 

m  (1  +  m/M)Z  mC2  +  2E 

b  =  2Z2e2/Mc2B2 

The  m,  v,  B,  y»  and  E  are  proton  parameters,  li  and  M  are  for 
Si28,  and  Em  is  the  maximum  energy  transfer  possible.  For  an 
energy  transfer,  Ep,  of  100  keV  or  greater  resulting  from 
100-MeV  proton  elastic  scatter,  oe  equals  208  mb.  This  is  a 
value  much  larger  than  the  p,  a  and  p  spall  cross  sections 
and  implies  that,  for  soft  error  thresholds  of  100  keV, 
elastic  scattering  will  contribute  to  soft  errors.  It  should 
be  noted  that  no  devices  currently  available  can  be  triggered 
by  so  little  energy,  but  100  keV  is  within  the  range  of 
projection  for  VLSI. 
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Summary 

1)  The  interaction  of  the  heavy  cosmic  rays  with  satel¬ 
lite-borne  microelectronics  can  be  viewed  along  two  lines: 

Soft  errors  (transient)  and  permanent  damage.  Soft  errors 
can  be  avoided  or  reduced  in  number  by  several  techniques; 

a)  design  in  accordance  with  Eq.  (7),  b)  fault  tolerant  cir¬ 
cuitry,  c)  redundant  programing.  No  consensus  on  the  event 
rate  at  which  soft  errors  become  a  "problem"  exists  at  this 
time.  Permanent  damage  in  general  cannot  be  shielded  against 
but  should  not  become  a  serious  problem  until  device  sizes 
reach  VLSI.  Failure  levels  can  only  be  approximated  until 
an  experimental  data  base  can  be  obtained  from  actual  circuits. 

2)  Generation  of  soft  errors  from  cosmic  proton  pri¬ 
maries  by  ionization  loss  along  single  tracks  is  not  seen  as 
a  possibility  for  LSI/VLSI. 

3)  Proton  generation  of  soft  errors  via  elastic  scat¬ 
tering  reactions  in  Si 28  become  increasingly  evident  when 
protons  of  energy  of  a  few  tens  of  MeV  are  present  and 
especially  when  the  threshold  energy  for  soft  error  reaches 
-lOOKeV.  Special  notice  should  be  taken  of  the  vulnerability 
during  solar  flare  time. 

4)  The  number  of  nuclear  reactions  which  will  occur  in 
volumes  as  minute  as  those  characteristic  of  VLSI  is  very 
small.  This  is  compensated  by  the  increased  number  of  ele¬ 
ments  per  chip  in  VLSI.  Thus,  upsets  caused  by  p,  a;  p  spall 
etc.,  reactions  can  still  be  anticipated.  If  the  population 
of  recoiling  residual  nuclei  which  result  from  inelastic 
scatter  and  which  have  energies  greater  than  a  few  hundred 
keV  is  large,  then  nuclear  inelastic  scatter  can  contribute 
an  event  rate  of  concern.  (During  the  publishing  time  of 
this  article  experimental  verification  of  the  proton  initia¬ 
ted  upsets  has  occurred. 20  evidence  suggests  that  in¬ 
elastic  scatter  can  produce  events  at  LSI  dimensions  as  well 
as  at  VLSI.) 

5)  The  electronics  package  vulnerability  assessment 
will  require  an  analysis  based  on  proton  Inelastic  scatter 
products  coupled  with  determination  of  the  appropriate  track 
distributions.  Experimental  and  theoretical  studies  of 
these  processes  were  undertaken  in  the  past2l»22  when  semi¬ 
conductors  were  much  larger.  The  process  needs  to  be  repeat¬ 
ed  using  current  and  projected  LSI/VLSI  circuitry  to  assess 
susceptibility  to  solar  flare  or  trapped  belt  fluxes. 
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Conclusion 

The  response  of  VLSI  circuitry  to  the  particle  fluxes 
encountered  in  Earth  satellites  can  be  partially  anticipated 
from  experience  in  nuclear  survivability.  Specifically,  the 
knowledge  of  device  behavior  when  exposed  to  fluxes  of  weakly 
ionizing  particles  which  cause  charge  trapping  in  the  oxides 
of  MOS  devices  applies  here.  Soft  error  rates  can  be  antici¬ 
pated  based  on  present  experience  with  LSI  memories  in  satel¬ 
lites  or  exposed  to  accelerator  fluxes  in  the  laboratory  and 
on  analyses  of  the  kind  presented  here.  The  possibility  of 
surprises  exists  because  of  permanent  damage  by  heavy  ion 
cosmic  rays  and  the  dramatically  nonuniform  dose  distribution 
surrounding  those  ion  tracks.  Additional  experimental  work 
also  is  needed  in  the  area  of  proton-nuclear  reactions  as 
they  pertain  to  the  creation  of  highly  ionizing  reaction 
products  in  electronics  materials. 
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